Search for a strongly decaying neutral 
charmed pentaquark 



The FOCUS Collaboration 



J. M. Link a P. M. Yager a J. C. Anjos b I. Bediaga b 
C. Castromonte b A. A. Machado b J. Magnin b A. Massafferri b 
J. M. de Miranda b I. M. Pepe b E. Polycarpo b A. C. dos Reis b 
S. Carrillo c E. Casimiro c E. Cuautle c A. Sanchez-Hernandez 
C. Uribe F. Vazquez L. Agostino d L. Cinquini d 
J. P. Cumalat d B. O'Reilly d I. Segoni d K. Stenson d 
c^, | J. N. Butler 6 H. W. K. Cheung 6 G. Chiodini I. Gaines 

^ ; P. H. Garbincius L. A. Garren E. Gottschalk P. H. Kasper 6 

A. E. Kreymer e R. Kutschke M. Wang L. Benussi f 
M. Bertani f S. Bianco f F. L. Fabbri f S. Pacetti f A. Zallo f 
M. Reyes g C. Cawlfield h D. Y. Kim h A. Rahimi h J. Wiss h 

R. Gardner 1 A. Kryemadhi 1 Y. S. Chung j J. S. KangJ 

B. R. Ko j J. W. Kwak j K. B. Lee j K. Cho k H. Park k 
G. Alimonti £ S. Barberis £ M. Boschini^ A. Cerutti^ 

P. D'Angek/ M. DiCorato^ P. Dmi e L. Edera £ S. Erba £ 
; P. Inzani^ F. Leveraro^ S. Malvezzi^ D. Menasce £ 

M. Mezzadri^ L. Moroni^ D. Pedrini^ C. Pontoglk/ F. Prelz^ 
M. Rovere^ S. Sala^ T. F. Davenport III m V. Arena 11 G. Boca 11 

G. Bonomi n G. Gianini 11 G. Liguori 11 D. Lopes Pegna 11 
M. M. Merlo 11 D. Pantea 11 S. P. Ratti 11 C. Riccardi 11 P. Vitulo 11 
C. Gobel H. Hernandez p A. M. Lopez p H. Mendez p A. Paris p 

J. Quinones p J. E. Ramirez p Y. Zhang p J. R. Wilson q 
T. Handler 1 " R. Mitchell 1 D. Engh s M. Hosack 8 W. E. Johns 8 
E. Luiggi 8 J. E. Moore 8 M. Nehring 8 P. D. Sheldon 8 
E. W. Vaandering 8 M. Webster 8 M. Sheaff* 

a University of California, Davis, CA 95616 
h Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, RJ, Brazil 
C CINVESTAV, 07000 Mexico City, DF, Mexico 



Preprint submitted to Elsevier Science 



7 February 2008 



d University of Colorado, Boulder, CO 80309 
e Fermi National Accelerator Laboratory, Batavia, IL 60510 
f Laboratori Nazionali di Frascati dell'INFN, Frascati, Italy 1-00044 
g University of Guanajuato, 37150 Leon, Guanajuato, Mexico 
^University of Illinois, Urbana- Champaign, IL 61801 
^Indiana University, Bloomington, IN 47405 
J Korea University, Seoul, Korea 136-701 
^Kyungpook National University, Taegu, Korea 702-701 
^INFN and University of Milano, Milano, Italy 
m University of North Carolina, Asheville, NC 28804 
n Dipartimento di Fisica Nucleare e Teorica and INFN, Pavia, Italy 
° Pontificia Universidade Catolica, Rio de Janeiro, RJ, Brazil 
p University of Puerto Rico, Mayaguez, PR 00681 
q University of South Carolina, Columbia, SC 29208 
r University of Tennessee, Knoxville, TN 37996 
s Vanderbilt University, Nashville, TN 37235 
t University of Wisconsin, Madison, WI 53706 



Abstract 

We present a search for a charmed pentaquark decaying strongly to D^~p. Finding 
no evidence for such a state, we set limits on the cross section times branching 
ratio relative to D*~ and D~ under particular assumptions about the production 
mechanism. 
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1 Introduction 

At the dawn of the QCD era, Jaffe proposed the existence of bound (mass 
below threshold for strong decay) multiquark states including QQqq states 
and the H dihyperon [1] based on calculations using the bag model [2]. Cal- 
culations a decade later [3,4,5] indicated that a S — — 1 charmed pentaquark 
(csqqq) is more likely to be bound than either the H dihyperon or a S = 
charmed pentaquark (cqqqq) due to the lack of a quark exchange process in 

1 See http://www-focus.fnal.gov/authors.html for additional author information. 
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the lowest decay channel, D~N 2 . Additional calculations using bag models 
and group theory assuming a one-gluon exchange interaction suggested that 
S — 0, — 1 charmed pentaquarks were probably unbound due to SU(3) sym- 
metry breaking and the finite mass of the charm quark [6,7,8]. Calculations 
using an instanton model indicated a nearly bound S — — 1 charmed pen- 
taquark [9] while calculations based on Goldstone boson exchange predicted 
unbounded S — — 1 charm pentaquarks but allowed bounded or unbounded 
S = charm pentaquarks [10,11]. Finally, the Skyrme (chiral soliton) models 
predicted deeply bound charmed pentaquark states with S — 0, — 1 [12,13,14]. 
Searches in 1998 and 1999 found no evidence for a weakly decaying charmed 
S — — 1 pentaquark from tt~N interactions [15]. 

The pentaquark field underwent a transformation between January 2003 and 
March 2004, when no less than ten independent pentaquark observations at 
a mass around 1540 MeV/c 2 were reported [16], with a presumed quark con- 
tent of (suudd). This triggered a rash of theoretical explanations, including 
predictions of charmed pentaquark states based on various quark models [17], 
Skyrme models [18], and lattice QCD [19] which predicted ®° c (cuudd) masses 
of 2600-3000 MeV/c 2 , 2700 MeV/c 2 , and 3450 MeV/c 2 , respectively, compared 
to the DN threshold of ~ 2805 MeV/c 2 . In March 2004, the HI Collaboration 
reported evidence for a S = charmed pentaquark state decaying to D*~p 
at a mass of 3099 ± 3 ± 5 MeV/c 2 and a statistical significance of 5.4 a or 
6.2 cr, depending on how the significance is calculated [20]. This precipitated 
additional theoretical effort into pentaquark states with a heavy quark [21]. 

This letter presents a search for the Q° c {cuudd) pentaquark candidate. We 
search using the same decay mode as HI (0°^ D*~p) and the other obvious 
decay mode, 0° — > D~p. Since the D* + statistics and data quality of the 
FOCUS experiment are much better than the observing experiment, and the 
production mechanism is similar, we should be able to confirm or refute the 
existence of the purported state. 



2 Event reconstruction and selection 

The FOCUS experiment recorded data during the 1996-7 fixed-target run at 
Fermilab. A photon beam obtained from bremsstrahlung of 300 GeV elec- 
trons and positrons impinged on a set of BeO targets. Four sets of silicon 
strip detectors, each with three views, were located downstream of the tar- 
gets for vertexing and track finding. For most of the run, two pairs of silicon 
strips were also interleaved with the target segments for more precise ver- 
texing [22]. Charged particles were tracked and momentum analyzed as they 

2 Charged conjugate states are implied unless explicitly stated otherwise 
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passed through one or two dipole magnets and three to five sets of multi- 
wire proportional chambers with four views each. Three multicell threshold 
Cerenkov counters, two electromagnetic calorimeters, and two muon detectors 
provided particle identification. A hadronic trigger passed 6 billion events for 
reconstruction. The average photon energy of reconstructed charm events is 
175 GeV. 

A candidate driven vertexing algorithm is used to reconstruct charm. In the 
case of D° — > K~ir + , K~it + it~it + (D + — > K~it + it + ), the charged tracks are 
required to verticize with CL > 2% (1%) with the correct summed charge. 
The momentum and vertex location are used as a "seed" track to find the 
production vertex, which must have CL > 1%. The D* + candidate is obtained 
from the decay D* + ^D°tt+. The soft pion must be consistent with originating 
from the production vertex and the track is refit using the production vertex 
as an extra constraint. 

The selection criteria were chosen to maximize S/ \fB. The signal, S, comes 
from a Monte Carlo simulation. Since the production characteristics of the 
pentaquark are unknown, we use PYTHIA to generate charm events with a 
S°. The S° is produced at a mass of 2.47 GeV/c 2 but is promoted to a mass 
of 3.1 GeV/c 2 before decaying to the decay mode of interest (with zero life- 
time). The background, B, is obtained from the wrong sign data (£>W+p) over 
the entire mass range of study (threshold to 4 GeV/c 2 ). This is an unbiased 
method for determining the selection criteria. 

Separating charm from hadronic background is primarily accomplished by re- 
quiring the decay vertex be distinct from the production vertex. A cut of 
II a t > 2, 4, 6 is applied to D° -> K~7r + , D° -> 7^+7^+, D+ -> K-it + it+ 
where £ is the distance between the two vertices and ag is the calculated 
uncertainty ((oy) ~ 500 fxm). Since hadronic reinteractions can fake a de- 
cay, requiring the decay vertex be located outside of target material reduces 
background. The out-of-material significance <r ut is positive (negative) for a 
vertex outside (inside) material. We require 2£/ae + max (—2, (Tout) > 4, 8, 12 
for D° -> K~ir + , D° -> 7^+7^+, D+ -> R-ti + ti + . To ensure the D° de- 
cay tracks do not originate from the production vertex, a cut is made on 
the change in production vertex confidence level when any D° decay track 
is added to the vertex (ACL pri ). We require ACL pri < 30% (0) for D° (£>+) 
decays. Also, the transverse momentum of the D candidate with respect to 
the line-of-flight obtained from the vertices must be less than 0.7, 0.4, and 
0.3 GeV/c for D° -> K'n + , D° -> K~it + it-it + , and D + -> R-ti + ti + , respec- 
tively. Since the signal contains a charged kaon, information from the three 
Cerenkov counters effectively suppresses backgrounds. The Cerenkov identifi- 
cation algorithm [23] returns negative log-likelihood (times two) values Wi(j) 
for track j and hypothesis i G {e, 7r, K, p}. In practice, differences in log- 
likelihoods between hypotheses are used such as A WW = W w — Wk- For the 
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Fig. 1. The normalized mass plots of D° ^ K~it + and D° — > K^tt + tt^tt + candidates 
(left) and D + — > iT _ 7r + 7r + candidates (right) are fit with a single Gaussian for the 
signal and a quadratic polynomial for the background. The energy release plot for 
D* + candidates of D° events (middle) is fit with a double Gaussian for the signal 
plus a threshold function aQ 1 ! 2 + (3Q 3 ^ 2 for the background. Events inside the 
vertical lines are selected for analysis. 

D*~ candidates, we require AW nK (K) > 0.5, W m i n (vr) - W n (n) > -5, and 
W min (7r s ) - W 7r (7r s ) > -5 where W min = min {W ie{e ^,K, P })- For the D + candi- 
dates, we require AW n K(K) > 2 and yV m i n (7r) — W 7r (7r) > —2. The remaining 
selection criteria are very efficient and mildly suppress some backgrounds. For 
the D° — > K~n + decay, the momentum asymmetry, , must be less 

than p ^ D 2oo^v/^ c ■ The summed "p\ of D daughters with respect to the D mo- 
mentum vector must be greater than 0.15 GeV 2 /c 2 . These last two cuts favor 
a decay of a heavy particle over combinatoric background. The charm signals 
obtained with these cuts are shown in Fig. 1. We select D° and D + candi- 
dates with an invariant mass within 2.5 <jm of the nominal mass where om is 
the calculated error on the mass. We select D* + candidates within 2 MeV/c 2 
of the nominal energy release, Q = M(D* + ) — M(D°) — m n +. These values 
were also chosen in the same optimization procedure. The proton candidate 
must be consistent with originating from the production vertex and must be 
strongly favored to be a proton by the Cerenkov systems; AW np (p) > 6 and 
AWk p (p) > 2. Again, these criteria were selected in the unbiased manner 
described above. 



3 Pentaquark search results 



The D*~p and D~p invariant mass plots, separated by charge states and 
summed together, are shown in Fig. 2 with fitted background curves superim- 
posed. The individual charge states show approximately an equal number of 
events in each state and so for the remainder of the analysis we combine the 
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Fig. 2. Invariant mass plots of D*~p candidates (left) and D~p candidates (right). 
The top and middle plots show the Dp and Dp contributions, respectively. The 
bottom plot includes both charge conjugate states. Plots are fit to a background 
shape: aQ b exp (cQ) where Q = M(Q C ) — M(D^~) — m p is the energy release. 

charge conjugate states. No evidence for a pentaquark at 3.1 GeV/c 2 or at any 
mass less than 4 GeV/c 2 is observed. To set a limit on the yield we first need in- 
formation about the width of the state. HI reports a measured Gaussian width 
of a T = 12 ± 3 MeV/c 2 with an expected resolution of a R = 7 ± 2 MeV/c 2 . 

We define the "natural" width as = \]&t ~ a \- Since <jt and <7r agree at 
better than 95% CL, the lower limit on a n is 0. Using the above measurements 
to construct a x 2 , we set x 2 — 3.84 to find the 95% CL upper limit on a^- 
The maximum value of for \ 2 = 3-84 occurs at = 16.6 MeV/c 2 with 
a R = 6.1 MeV/c 2 and <jt = 17.7 MeV/c 2 . We construct limits on yields under 
the two extreme assumptions about the pentaquark "natural" width, = 
and o"v = 16.6 MeV/c 2 . To a very good approximation, the upper limit in- 
creases monotonically with <jn so this provides the extreme range. In both 
cases, <ttv is added in quadrature to the experimental resolution to obtain the 
total width, ot- The experimental resolution increases linearly from 2 MeV/c 2 
at threshold to 14 (13) MeV/c 2 at 3.95 GeV/c 2 for D*~p (D~p). 

The D~p (D*~p) mass plots are fit in 1 MeV/c 2 steps from 2.84 (2.98) to 
3.9 GeV/c 2 . At each step, the mass and width of the Gaussian signal are fixed 
while the signal yield and background parameters are fitted. A binned log- 
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Fig. 3. Charmed pentaquark yields and upper limits. Left (right) plots show results 
for D*~p (O c — > D~p). Top (bottom) plots show results for a natural width of 
(16.6 MeV/c 2 ). The shaded region includes the 1 a errors with the central value 
in the middle. The outer curves show the upper and lower limits. 

likelihood fit using Minuit [24] is performed with a background shape given 
by aQ b exp (cQ) where Q is the energy release. The 1 a error is denned as the 
point where Alog£ = 0.50 relative to the maximum log£, while continually 
adjusting the background parameters to maximize log/}. The 95% CL lower 
limit is denned similarly with Alog£ = 1.92. Both are obtained using Mi- 
nos [24]. The 95% CL upper limit is constructed as follows. The likelihood 
function £ versus yield is determined by maximizing log£ for many different 
(fixed) yields, allowing background parameters to float. The likelihood func- 
tion is integrated from a yield of to oo. The 95% CL upper limit on the yield 
is defined as the point where 95% of the total likelihood is between a yield 
of and the upper limit. The fitted yield, l-a errors, and 95% CL limits are 
shown in Fig. 3. 

To compare with other experiments, the limits on yield must be converted to 
limits on production times (unknown) branching ratio. We choose to nor- 
malize the 9° production cross section to the D meson production cross 
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section from which it is reconstructed. That is, we attempt to determine 
a (Q° c ) ■ BR(Q° C ^ D-p) /a (D-) and a (0°) • BR (Q° c ^D*-p) /a (£>-). 

The FOCUS detector is a forward spectrometer and therefore acceptance de- 
pends on the produced particle momentum. The production characteristics of 
the charmed pentaquark is by far the largest source of systematic uncertainty 
in this analysis. We choose a particular production model to obtain limits and 
provide sufficient information about the experiment for other interested parties 
to obtain limits based on other production models. The production simulation 
begins with a library of e~ and e + tracks obtained from a TURTLE simula- 
tion [25] of the Wideband beam line. From this library, an individual track 
is drawn and bremsstrahlung photons created by passage through a 20% Xo 
lead radiator. The charm cross section is then applied to the bremsstrahlung 
photons using Monte Carlo rejection. A photon which produces a charm event 
is sent to the Pythia event generator [26]. 3 Since Pythia does not pro- 
duce charmed pentaquarks, another particle must be chosen to represent the 
charmed pentaquark. Charmed baryons are the natural choices since they have 
a single charm quark like the charmed pentaquark and are closer to the correct 
number of total quarks than charmed mesons. Since it is possible to adjust the 
mass of the chosen charmed baryon in Pythia, it is not necessary to pick the 
charmed baryon with the highest mass. Other than mass, the most important 
effect on the production is the number of quarks a particle has in common 
with the initially interacting hadrons, due to the nature of the Pythia string 
fragmentation model. The S° and X+ particles are chosen to represent the 
extremes in the production of a charmed pentaquark. Other than the charm 
quark, the E^(csd) (E+(cuef)) can obtain at most 50% (100%) of the remaining 
quarks from the target nucleon valence quarks, while the Q° c {cuudd) can take 
75%. In all cases, the charge conjugate particles must obtain all quarks from 
the vacuum. The mass of the particle chosen to represent the pentaquark, S° 
or E+, is set to the appropriate value in Pythia, by setting PMAS(190,1) or 
PMAS(187,1), respectively. This method differs from that used to generate the 
Monte Carlo sample for cut optimization since in this case the produced mass 
is changed to the appropriate value. For a photon which passes the charm 
cross section, the Pythia generator is run up to 1000 times searching for 
a S° This fails approximately 90% of the time, especially for low en- 

ergy photons. When this happens, a new charm producing photon is selected. 
This changes the photon spectrum for pentaquark producing events. The ini- 



3 Pythia version 6.127 is used with modifications to many parameters which pro- 
vides a better match to the FOCUS charm data: MSTP(14): 20^0, PARP(2): 10 
3, PMAS(4,1): 1.35 -► 1.6, MSTJ(21): 2^0, PARJ(l): 0.1 -► 0.3, PARJ(3): 0.4 -► 
0.25, MSTJ(12):2 -► 4, MSTP(92):4 -► 2, PARP(96):3 -► 2, PARP(97):1 -► 2, 
PARP(91): 0.44 ->0.6, PARP(93): 2 -► 3, PARP(99): 0.44 -► 0.6, PARP(100): 2 -► 3, 
PARJ(21): 0.36 -^0.425, MSTJ(13):0-»1, PARJ(41): 0.3 -► 0.25, PARJ(42): 0.58 -► 
0.7, PARJ(45): 0.5^0.75, PARJ(36): 2->l. 
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Fig. 4. The initial bremsstrahlung spectrum peaks at low energies. Applying the 
rising charm cross section results in a harder photon spectrum for charm inducing 
photons. The photons producing 3.1 GeV/c 2 H c and S c particles are harder than 
generic charm inducing photons except for a low energy peak for S c production due 
to the charm quark combining with quarks from the target nucleon. 

tial bremsstrahlung spectrum, applied charm cross section, charm producing 
photon spectrum, and pentaquark producing photon spectra for S° and 
are shown in Fig. 4. After production, the H° or £+ is changed to the 0° 
with zero lifetime and forced to decay in the mode of interest. The momen- 
tum of the generated particles, the momentum of the reconstructed particles, 
and the ratio of the two (efficiency versus momentum) are shown in Fig. 5 
for a pentaquark mass of 3.1 GeV/c 2 . Since there are many low momentum 
pentaquarks from E+ which are not reconstructed, the overall efficiency for 
pentaquarks produced as is lower than for pentaquarks produced as H°. 
The low generated momentum spike for pentaquarks produced as E+ is a re- 
sult of the produced charm quark combining with a u and a d quark from the 
target. The efficiency versus mass is shown in Fig. 6 where the top (bottom) 
curve shows the efficiency for pentaquarks produced as 5° (£+). 

The limits on 0° production are made relative to £)M~ production. Therefore, 
the efficiency for D*~ and D~ must be determined. Pythia is again used to 
model charm production. In this case we have verified that Pythia correctly 
models charm production in the region over which we have acceptance. This 
is less than 10% of the total charm produced. Note, however, that only the 
ratio of efficiencies between 0° and production need be determined. 

Therefore, discrepancies between Pythia and reality in the low momentum 
(unobserved) region are likely to cancel. The global efficiency for D*~ (D~) 
reconstruction is 2.81% (3.05%). 
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Fig. 5. All four plots show generated and reconstructed events (histograms) and the 
calculated efficiency (points with error bars) versus pentaquark momentum. The 
top (bottom) plots are for pentaquarks produced as H° (£+). The left (right) plots 
are for the decay D*~p (Q®^D~p). 
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Fig. 6. Efficiency (including acceptance) versus mass for reconstructing @®^>D*~p 
(left) and 0j? — > (right). The top (bottom) curves are obtained from pen- 

taquarks produced as (£+) while the middle curve is the average. 



We attempt to determine: 

a(Q° c ^Dp) -BR (&° c ^ Dp) Y e o 



a(D) 



Y D 



10 



where D is either D*~ or D~ . The D yield (Yd) and pentaquark yield {Y@o) are 
obtained from Fig. 1 and Fig. 3, respectively. The uncertainty on the D yield 
is less than 1%, insignificant compared to other uncertainties in the analysis. 
The remaining quantity is the ratio of efficiencies, e£>/e o. Our uncertainty in 
€d is much less than eeo since the Monte Carlo simulation uses a tuned version 
of Pythia which accurately reproduces the D production for the region in 
which the FOCUS detector has acceptance. The dominant source of system- 
atic uncertainty comes from a lack of knowledge of the pentaquark production 
characteristics which is reflected in e o. An estimate of this systematic un- 
certainty is obtained from Fig. 6 which shows the efficiency for 0° produced 
as S° (top curve), (bottom curve), and the average (middle curve). We 
take the average as our central value for the efficiency (eo) and the difference 
between the average and top (or bottom) curve as the 1 a error on the effi- 
ciency (cr tQ ). To include the error we fit directly for the corrected yield Y e o/eQO, 
include the efficiency as a free parameter and use a modified log-likelihood: 
log/y = log£ — 0.5 (efit — e ) 2 /of where 6fi t is the fitted efficiency which is al- 
ways equal to eo when log C is maximized. The corrected yield, 1 a errors, and 
95% CL upper limit is obtained in the same manner as described previously 
for the uncorrected yield and is shown in Fig. 7. 



4 Conclusions 



FOCUS has published the most precise measurements of the mass and width 
of four D** states [27] with decay modes similar to that of the charmed pen- 
taquark. However, from a FOCUS data set of 104,000 D*~ and 137,000 D~ de- 
cays, we find no evidence for a charm pentaquark decaying strongly to D^~p 
over the mass range of 30 MeV/c 2 above threshold to 3.9 GeV/c 2 . Results for 
yields and cross sections versus mass are found in Figs. 3 and 7. A summary 
of these limits over the entire mass range is given in Table 1. The sample of 
D*~ events is more than 30 times larger as well as cleaner than the sample 
used by HI in the paper which found evidence for the charm pentaquark [20]. 
While HI finds -1% of the D*~ from 0°, FOCUS sets a limit of < 0.075% at 
95% CL. In addition, the production is similar between the two experiments; 
virtual (real) photons on protons (nucleons) for HI (FOCUS). Thus, the HI 
result is either a statistical fluctuation or the result of an unusual production 
mechanism which increases the charm pentaquark to charm cross section by 
a factor of at least 10 in HI relative to FOCUS. A simple study to determine 
the extent of the difference was performed. Pythia was used to generate D* + 
and O°(3100) events using production appropriate for FOCUS and HI. The 
O°(3100) was modeled as the average of S° and E+, both at 3.1 GeV/c 2 . Only 
events in the acceptance of the experiment were accepted by requiring the pen- 
taquark momentum be greater than 50 GeV/c for FOCUS and requiring HI 
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Fig. 7. Cross section times branching ratio results relative to single D production. 
Left (right) plots show results for C ^ D*~p (B c ^ D~p). Top (bottom) plots show 
results for a natural width of (16.6 MeV/c 2 ). The shaded region encompasses the 
la errors with the central value in the middle. The top curve shows the upper limit. 
Both include systematic uncertainties. 



events to have 1<Q 2 <100 GeV 2 /c 2 and 0.05 <y< 0.70. From these samples, 
the production rate of O°(3100) to D* + was found to be 1.58% (1.19%) for 
HI (FOCUS). Scaling by the observed production rate of 1% at HI, the pro- 
duction rate of O°(3100) relative to D* + should be (1.19/1.58) ■ 0.01 = 0.75% 
for FOCUS. For this FOCUS analysis, the O°(3100) efficiency relative to D*+ 
is 54%. Therefore, from the FOCUS sample of 107,525 D* + events, one would 
expect 107525 • 0.0075 • 0.54 = 435.5 O°(3100) events. This signal has been 
superimposed on the FOCUS data in Fig 8. This signal would be easily ob- 
servable. The non observation by FOCUS is also consistent with published 
results from ALEPH [28] and ZEUS [29]. 



12 



Table 1 

Summary of upper limits (UL) on pentaquark yields and cross sections including 
systematic uncertainties. Results represent the maximum UL over the mass range 
shown in Figs. 3 and 7. 
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Fig. 8. The FOCUS data (from the left plot of Fig 2) plus a simulated Gaus- 
sian signal of 435 events based on the observed production at HI. The top plot 
shows the signal with a width due to experimental resolution (a = 4.15 MeV/c 2 ) 
only while the bottom plot shows the signal with a width due to the experimental 
resolution plus the maximum natural width (at 95% CL) based on the HI data 
(a = 4.15 16.6 = 17.1 MeV/c 2 ). 
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